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The dynamical behaviour of many-body systems is often richer than what can be anticipated from 
their static properties'^. Here we show that in closed quantum systems this becomes evident by 
considering time-integrated observables as order parameters. In particular, the analytic properties of 
their generating functions, as estimated by full-counting statistic s 16 " 19 *, allow to identify dynamical 
phases, i.e. phases with specific fluctuation properties of time-integrated observables, and to locate 
the transitions between these phases. We discuss in detail the case of the quantum Ising chain in a 
transverse field 2 . We show that this model displays a continuum of quantum dynamical transitions, 
of which the static transition is just an end point. These singularities are not a consequence of 
particular choices of initial conditions or other external non-equilibrium protocols such as quenches 
in coupling const ant They can be probed generically through quantum jumpP^' statistics of 
an associated open problem, and for the case of the quantum Ising chain we out line a possible 
experimental realisation of this scheme by digital quantum simulation with cold ion s ' 22 * 23 ! 



In order to identify and classify phase transitions in a 
many-body system it is necessary to determine the ap- 
propriate order parameters, i.e. system extensive observ- 
ables which allow to distinguish between phases, and the 
corresponding conjugate fields which can drive the sys- 
tem across the transitiorP^. In the case of dynamical 
phase transitions it is often assumecP^ that it suffices to 
study properties of the steady state and that stationary 
observables can function as faithful witnesses of drastic 
changes in the dynamics. While it is certainly true that 
if one drives a system across a static phase transition it 
is likely that dynamics will undergo a transition at the 
same pointP^ the converse may not be true in general, 
as dynamical transitions can also occur away from static 
one!™. 

In many-body systems, therefore, dynamics cannot al- 
ways be inferred from statics. To uncover the full range 
of dynamical phase behaviour it is necessary to consider 
strictly dynamical observables. In the case of both classi- 
cal and open quantum systems, thermodynamic or static 
transitions relate to singular changes in ensembles of 
configurations or stated^, while dynamical o nes re late 
to singular changes in ensembles of trajectories^^. In 
classical systems these trajectories are histories of time 
evolution in configuration spaced for open quantum 
systems they correspond to the tim e record of quanta 
emitted by the system into the batr l 12 | 2Q ^4 Appropri- 
ate order parameters are time-integrated observables as 
they capture the dynamical fluctuations in these trajec- 
tories which give rise to dynamical transitions. Dynam- 
ics can t hen b e studied by considering the full counting 
statistic j 12 * 16 ^ ^ 2 ^ of such time-integrated observables. 
In the long-time limit this FCS approach, supplemented 
with large-deviation method^, yields quantities akin to 
free-energies for ensembles of trajectories, whose analytic 
properties in terms of "counting" fields (the fields mathe- 
matically conjugate to dynamical obser vables) rev eal dy- 
namical phase structure and transition s 1 1 l 12 l 18 iMJ j 

In the real time dynamics of closed quantum systems 



we do not have a similar concept of observable trajec- 
tories. Nevertheless, we show here it is possible to pur- 
sue a strategy analogous to that of open problems by 
studying the properties of generating functions of mo- 
ments or cumulants of time-integrated quantities, objects 
which are well defined in closed quantum systems. Sin- 
gularities of these generating functions in the long-time 
limit locate quantum dynamical phase transitions. Fur- 
thermore, these can occur for values of static parameters 
away from those of static transitions. As a consequence, 
even far from static transitions, fluctuations associated 
to dynamic singularities will become manifest in time- 
correlation functions and therefore directly influence the 
observed dynamics. 

Consider a closed quantum system with (time- 
independent) Hamiltonian H. A generic dynamical order 
parameter is given by the time-integrated observable, 

K t = f dt'k{t'), (1) 
Jo 

where k(t) = U}kU t is an operator in the Heisenberg pic- 
ture, Ut = e~ ltH the evolution operator (we set h = 1), 
and fct = k. Dynamical fluctuations are captured by the 
expectation value of K t and its higher order cumulants. 
To obtain these we define the moment generating func- 
tion (MGF) of K t : 

Z t (s) = (Tl(s)T t (s)), (2) 
where the modified evolution operator is defined as 

T t (s) = e- it6 -, H s = H-^k (3) 

It is easy to see from these definitions that indeed Z t (s) 
generates the moments of K t through its derivatives, 
(K?) = (-) n d s Z t (s)\ s ^ , while the logarithm of the 
MGF, 6t(s) = log Z t (s), is the cumulant generating func- 
tion (CGF). The expectation value denoted by (•) could 
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be over either a pure state or a mixed state. If a sys- 
tem has dynamical phase transitions they become man- 
ifest in singularities of these generating functions. The 
defini tions p| )-([3|) are a form of full counting statistics 
(FCS) 16 | 17 | 19 | for K t . In contrast to the standard FCS ap- 
proach we consider s real. This makes H s non-Hermitian 
and Tt(s) non-unitary, which in turn implies that 



0(s) = lim -G t (s) 

t— >oo t 



(4) 



is well defined. This am ounts to extending the traject ory 
method of RefsP^ 2 * 18 * 24 * (sometimes called a s-en sembl e"^) 
to closed quantum systems. In open systend 2 ^!!, 6(s) 
corresponds to the large-deviation 2 -^ rate function of the 
MGF, that is, Z° pen (s) « e te ^ for t large, and there- 
fore plays the role of a free-energy for ensembles of 
trajectorie d 12 * 24 *. In closed quantum systems such prob- 
abilistic interpretation not always possible 27 . But just 
like in open problems, 0(s) determines the dynamical 
phase structure, and specifically its singularities as a 
function of s locate quantum dynamical phase transi- 
tions where the cumulants of K t change in a singular 
way. The function Q(s) is the key quantity to compute in 
this approach. We now illustrate these ideas by study- 
ing dynamical transitions in time-integrated observables 
of the one-dimensional quantum Ising model, and show 
how such transitions can be probed through quantum 
jump statistics of an associated open problem. 

The quantum Ising model in a transverse field in one 
dimension has a Hamiltonian 2 , 
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We assume periodic boundary conditions, and define en- 
ergy in units of the exchange interaction between spins. 
This is the prototypical model for a system displaying a 
quantum phase transition 2 , which occurs in the limit of 
N — >• oo at |A C | = 1, from a disordered state for |A| > 1 
to an ordered state for |A| < 1. 

We choose as a time-integrated observable, Eq. ([!]), 
the time-integral of the transverse magnetisation, k = 
af, and we study for convenience the case when 
the system is in its ground state |0). The MGF 
of Eq. (|]) is Z t (s) = (0\T}(s)T t (s)\0), where the 
non-Hermitian Hamiltonian (K3| is given by H S (X) — 
H(X + is/2). Thus defined, Z t (s) can be calcu- 
lated exactly using standard free-fermion techniques 2 
(see Methods). The eigenvalues of H s are complex, 



€k(s) = 2y (A — cos k + is/2) 2 + sin 2 fc, which guaran- 
tees the time-convergence in Eq. Q. From their imagi- 
nary parts we get 6(s) = 2 ^2 k>0 \Tm €k(s)\. In the large 
size limit the corresponding per-spin function, 0(s) ee 
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FIG. 1. Phase transitions in time- integrated observ- 
ables of the quantum Ising model, a, Dynamical phase 
diagram: the dynamically disordered and dynamically or- 
dered phases are separated by a curve of second-order transi- 
tions corresponding to the singularities of of 9(s) on the semi- 
circle dT>; the black circles on the A axis indicate the location 
of static transitions; the regions I, II, III are defined by the 
structure of the state \s) (see main text and Methods); yellow 
dashed lines indicate the cuts plotted in the next two panels, 
b, Dependence on s for fixed A of the average time-integrated 
transverse magnetisation k, s (green) and its dynamic suscep- 
tibility Xs (black) which diverges as the phase boundary is 
approached from inside V; also shown is the s-biased static 
magnetisation m xs (green) which is directly related to 0(s), 
see main text, c, The same as in (b) but as a function of A 
for fixed s. 



limTv^oo^V 1 0(s), becomes, 



6(8) = { 



ilm{(l + X s )[E(^) 



where E(x) and E((p, x) are the complete and incomplete 
elliptic integrals of the second kind, A s ee A + is/ 2, and 
k\ = cos -1 A. The region V in the (A, s) plane is that in- 
terior to dV = {(A, s) : A 2 + (s/2) 2 = 1}. The expression 
above is for s > as 0(s) = 0(—s). 

The function 0(s) encodes properties of the long-time 
behaviour of the MGF Z t (s) and through it that of the 
cumulants of the time-integrated observable K t . Specif- 
ically, its analytic properties, which are given by the re- 
sponse of the spectrum of H to the deformation of Eq. 
([3]), determine the possible dynamical regimes, or dynam- 
ical phases, of the system. The corresponding phase di- 
agram is shown in Fig. 1(a). For each value of A there 
are two dynamical phases separated by a phase transi- 
tion at s c (X) ee 2sink\. The transition is second order 
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FIG. 2. Dynamical transitions of the closed Ising model from auxiliary open model via digital simulation with 
cold ions, a, World lines of 7 ions simulating a 6-spin Ising ring with dissipation. Shown are the gate operations which 
evolve the system over a single time step. Single-ion operations (blue squares) corresponding to interaction with the transverse 
field and two-ion operations (orange squares) simulating the exchange interaction determine the coherent evolution. Two-ion 
gates with an ancilla ion (green squares) simulate the dissipative dynamics. Measurement of the ancilla ion (red diamonds) 
determines whether a quantum jump occurs (see Methods), b, Theory and numerical simulation of Po(t) for a 6-spin Ising 
ring in its ground state with A = cos(57r/6) for (dimensionless) time steps of length 0.1 and 0.01 (labelled). The equivalent 
simulations with spins initialised in the state |— ) are shown with the label (*). c, Simulations of Po as a function of A and s/2 for 
4, 6 and 8 ions at t = 5, 5 and 8 respectively. Green overlays show the unit circle with special positions (s/2, A) = (sin k, cos k) 
(see main text). 



at all points on the curve s c (A). The expected value 
of the time-integrated observable K t (per unit time), 
n s = —6 f (s), is the order parameter of the dynami- 
cal phases. It is continuous at s c (A), while its corre- 
sponding susceptibility, Xs = diverges when s c (A) 
is approached from inside D\ see Figs. l(b,c). This is 
a square root singularity, \s ^ \s — s c (A)| -1 / 2 , except 
at (A, s) = (±1,0) where it is logarithmic. This lat- 
ter critical behaviour is the same of the static transverse 
susceptibilt^P: the static phase transitions are the end- 
points of the dynamical singularites. 

The nature of the dynamical phases can be understood 
from the state \s) = lim^oo T t (s)\i). For the case we are 
considering the initial state is the ground state, \i) = |0), 
but in the long-time limit \s) does not depend, up to 
a normalisation, on this initial condition (see Methods). 
Similarly, we define the s-biased expectation values of 
observables by (6) 8 = lim^ Z- 1 (s)(0\T^s)OT t (s)\0). 
From Eqs. Q-(|3]) it is easy to see that this expectation 
value taken for the operator k is directly related to the 
long-time CGF, (k) s = —0(s)/s. In our case k is the 
transverse magnetisation, so that m xs = TV -1 J2 i (crf ) s 
has the same singular behaviour as see Fig. l(b,c). 

For all values of (A, s) outside and inside the region X>, 
the states \s) are smoothly connected, but they change in 
a singular manner across the boundary dV. We denote 
these two regions dynamically disordered and dynami- 



cally ordered, respectively, since the spectra of H S (X) are 
smoothly connected to the spectra of H(\) which define 
the corresponding disordered and ordered static phases; 
see Fig. 1(a). The distinction between the dynamical 
phases is the following: if | A| > 1, the system is in the dis- 
ordered static phase, and there is no singular behaviour 
in its real-time dynamics; while if |A| < 1, corresponding 
to the ordered static phase, dynamical fluctuations will 
display singular behaviour, manifestly in the cumulants 
of time-integrated observables, due to the singularities 
of their generating functions at s c (A). It is important 
to note that the transitions our method reveals do not 
depend on non-equilibrium protocols^', such as quench- 
ing across a static phase boundary, or on a particular 
choice of initial state, but are an intrinsic feature of the 
spectrum of the problem. 

We now show how to probe the dynamical transi- 
tions of closed quantum systems described above from 
the quantum jump statistics of an associated open prob- 
lem. This connection is due to the fact that the 
MGF Z t (s) for the time-integrated observable K t can 
be obtained from the waiting time distribution between 
quantum jump^EU Q f an auxiliary open quantum sys- 
tem, since the H s evolves a density matrix p(t) ac- 
cording to p(t) = —i[H,p(t)] — §{fc,p(£)}. This is a 
Lindbla d 20 * 21 * master equation without recycling terms, 
to which we can associate an open quantum system 
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described by a full Markovian master equation p = 
-i[H,p] + Y,i ( L iP L l ~ H L l L ^P})' Specifically, we 
identify k with a set of quantum jump operators Li con- 
structed such that s ^ ji L\Li = sk. The operator T t (s) 
is then the same one that evolves the associated open 
quantum system between quantum jump events, with s 
being the decay rate of quantum jump processes. The 
MGF Z t (s) of the closed system then equals the proba- 
bility Po(t) that no quantum jumps occur up to a time 
t in the associated open system. This allows the MGF 
to be determined by preparing the closed system in an 
initial state and coupling it to an appropriate environ- 
ment. By finding the distribution of waiting times until 
the first quantum jump occurs, Po(t), the MGF of the 
close system can then be inferred. 

For the Ising model ([5| above, if we make a trivial shift 
such that k = \af + 1), we can choose quantum jump 
operators Li = y/2s \-)a(+\, where crf\±)i = ±\±)u 
and i runs over the sites of the lattice. Such open 
quantum system can be studied experimentally using the 
digital simulation techniques available in cold-ion sys- 
tem J 22 * 2 ^ . By applying a series of one- and two-ion gate 
operations, the open-system time evolution is approxi- 
mated by a Trotter decomposition of T t (s) with finite 
time steps, as sketched in Fig.|2|a). The implementation 
of dissipative dynamics involves the use of an ancilla ion 
whose state is measured after each dissipative gate oper- 
ation, Fig. [2|a), with the result determining whether a 
quantum jump has occured 28 . Repeating the experiment 
many times up to the first quantum jump allows Pq to be 
estimated; the MGF Z t (s) is extracted at different s by 
using different decay rates for the dissipative dynamics. 

We show simulations for Po(t) using the Trotter de- 
composition in Fig. |2^b), demonstrating that Po(t) can 
be found accurately at finite times. Figure [2^b) further 
shows that the behaviour of Po(t) at long times is rescaled 
to larger probabilities if all spins are initialised in the 
state |— ), rather than the Ising ground state, as this state 
is annihilated by the jump operators. This underlines the 
fact that Z t (s) encodes properties of the whole spectrum, 
not just the ground or low lying states, so that at long- 
times the precise nature of the initial state does not mat- 
ter. This allows features in the closed-system MGF to be 
explored at longer times with a smaller chance of each 
experimental run being terminated by the first quantum 
jump. Figure [2jc) shows the result of simulations us- 
ing this initial state for different decay rates s with dif- 
ferent magnetic fields A for N = 4, 6 and 8 ions. We 
see, even at finite times, marked features close to the 
semicircular transition line which exists in the thermo- 
dynamic limit. Each of these features lies close to posi- 
tions (A, s/2) = (cos fc, sin k) on the unit circle, where the 
values for k are the quasi-momenta associated with the 
excitation spectrum of the TV-spin Ising model. 

We have shown here, by explicitly extending the con- 
cept of an order parameter to the dynamical domain, 
that fluctuations in time-integrated observables reveal 



dynamical singularities in the quantum Ising model even 
away from its static transition, and that these can be 
probed in quantum jump statistics of an associated open 
problem. While these dynamical singularities are strictly 
present only in the limit of large size and time, we 
have shown that clear evidence of them can be observed 
in finite systems in regimes accessible to experiments. 
The approach we presented should help reveal dynam- 
ical quantum transitions that go beyond static ones in 
closed quantum many-body systems in general. Our 
work here should also connect to studies of thermalisa- 
tion in closed quantum system d 6 * 29 !, which often focus 
on time-integrated quantities under the assumption that 
they converge to expectation values of statistical ensem- 
bles: our results show that these quantities can fluctuate 
in a singular manner and this may strongly influence the 
ability of a system to thermalise. 

METHODS 

Diagonalisation of H s — The non-Hermit ian Hamil- 
tion is solvable via a Jordan- Wigner transformation 
and Bogoliubov rotation 2 which maps H s to a free 
fermion model with a complex dispersion relation, e k = 

2\J (A + y — cos k) 2 + sin 2 /c. These k modes are dis- 
crete and take values k = irn/N, where n = — N + 
1, — N + 3, . . . , N — 1, for the Ising chain with peri- 
odic boundary conditions. For specificity we focus on TV 
even. The diagonal form is H s = Yk e k( s ) (AkA^ — 1/2), 
where (Ak,Ak) are a conjugate fermion operator pair, 
{Ak,Ak f } = Sk,k', but because of the non-Hermitian na- 
ture of H s we have that A^ ^ A\. The key property for 
evaluating our Z (s,t) is that the vacuum of the unper- 
turbed Ising model, |0), may be expressed as a BCS state 
of the new non-Hermitian Hamiltonian. 

|0) = (g)[cosa||O fc ,0_,) s -isin^|l fc ,l_,)J (6) 

fc>0 

Here (g) stands for direct product, |n/ c ,n_/ c ) s indicate 
occupation states of the fermionic modes with |fc| that 
diagonalise H S1 and the coefficients are related to the 
Bogoliubov angles 2 , </>|, via a k = ^ k ~^ k ; where fa — 

0| =o . Requiring all off-diagonal terms in the H s to vanish 
we find these angles are given by tan fa k = sin k/[X + y — 
cos/c]. From this one may evaluate the partition sum 
directly, 

Z(*>t) = n(l cosa fe| 2cosh [ 2Xm K)] e " 2Xm(e " )£ ( 7 ) 

k>0 

+ | sinc^l 2 cosh [2 Imfaaie 21 ^'* 
+isinc4 cosa- s sinh [2 Xm(a|)]e- 2i ^ e(e ^ )£ 

-i sin a k s cos a s k sinh [2 Xm{a%)]e 2in < e8 ^ . 

Eq. ([7]) has a second-order singularity on the curve dV. 
Due to the connection between the quantum Ising chain 
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and the two-dimensional classical Ising modeP this criti- 
cal critical curve is related to the Lee- Yang zeros^ of the 
latter. 

The state |s), see main text, is written in terms of the 
fermionic modes |n/ c ,n_/ c ) s . The precise occupation of 
the fermionic levels depends on A. By applying T t (s) to 
the initial state ([6| we obtain in the long time limit, up 
to constants, 

r <g> fc >ol°^ -fc) s x > 1 

l*)ocj ® k >kj°k,0-k) a ® k<kx \lk,l-k) a -KA<1 

I ®; b >ol 1 *> 1 -fc> a A< - 1 

where for |A| < 1 the wavevector k \ is defined through 
A = cosk\. These are the regions (I, II, III) indicated in 
Fig. 1(a). We may now compute the s-biased expectation 
value of any observable. Specifically, the magnetisation 
m xs defined in the main text reads, 



1 

TV 
1 

N 



Xm(e s k )<0 



£ 

Xm(e*)>0 



V 



sin(f) 



cosh[22"ra(a|)] 



cos 



cosh[22ra(a|)] 



Digital Simulation — Digital simulation of open quan- 
tum systems is based upon performing Trotter decompo- 
sitions on the unitary operator generating coherent evo- 
lution, and introdu cing a n ancilla spin to simulate the 
dissipative dynamic d 22 * 28 *. The time evolution of the state 
of the system ancilla system, \ip) s <g)\A) a , is approximated 



by a series of unitary transformations corresponding to 
short time steps St. A single time step evolves the com- 
bined system-ancilla state according to 
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n^n 
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k=i 



N 

IT 



® |0)a- (8) 



The two rightmost terms evolve just the system accord- 
ing to its self Hamiltonian: single-body operations which 
result from the magnetic field A are applied before a se- 
ries of two-body operations describing the Ising spin-spin 
interaction; see Fig. 2(a). The gate operation Gj acts on 
the Hilbert space of the spin at site j and the ancilla spin, 
which is initially prepared in state |0) a , such that 

Gj\+)j <8> |0) s = cos^)|+) i <g> |0) a - zsin^l-^- <8> |l) a 

G j \-) j ®\0)s = \-) j ®\0) a 

with (j) = y/jSt. After each application of the gate G\ 
the state of the ancilla is measured in the |0) a , |l) a ba- 
sis, see Fig. 2(a) in the main article. The upper limit 
on the outmost product in Eq. ([8| is for n = N for all 
the time steps where it is applied, except for the final 
one where n<> < N. In the final time step the ancilla 
is measured in the state |l) a , indicating that a quantum 
jump process has occured; for our purposes, the simula- 
tion is now terminated and the time taken for a jump 
to occur is recorded. If the state |0) a is measured, the 
Trotterised evolution continues with the ancilla ion reset 
in state |0) a . For small </)« 1, the dissipative gate op- 
erations accurately simulate the evolution of a Lindblad 



master equation with jump operators Lj 



,<+i 



Acknowledgments. — The work was supported by EP- 
SRC Grant no. EP/I017828/1 and Leverhulme Trust 
grant no. F/00114/BG. 



1 Peliti, L. Statistical mechanics in a nutshell (Princeton 
University Press, 2011). 

2 Sachdev, S. Quantum Phase Transitions (Cambridge Uni- 
versity Press, 2011). 

3 Binder, K. & Kob, W. Glassy materials and disordered 
solids (World Scientific, 2011). 

4 Diehl, S., Tomadin, A., Micheli, A., Fazio, R. & Zoller, 
P. Dynamical phase transitions and instabilities in open 
atomic many-body systems. Phys. Rev. Lett. 105, 15702 
(2010). 

5 Henkel, M. & Pleimling, M. N on- Equilibrium Phase Tran- 
sitions: Volume 2 (Springer, 2010). 

6 Polkovnikov, A., Sengupta, K., Silva, A. & Vengalattore, 
M. Nonequilibrium dynamics of closed interacting quan- 
tum systems. Rev. Mod. Phys. 83, 863 (2011). 

7 Calabrese, P., Essler, F. H. L. & Fagotti, M. Quantum 
quench in the transverse-field ising chain. Phys. Rev. Lett. 
106, 227203 (2011). 

8 Dorner, R., Goold, J., Cormick, C, Paternostro, M. & Ve- 
dral, V. Emergent thermodynamics in a quenched quan- 



tum many-body system. Phys. Rev. Lett. 109, 160601 
(2012). 

9 Hinrichsen, H. Non-equilibrium critical phenomena and 
phase transitions into absorbing states. Adv. Phys. 49, 
815-958 (2000). 

10 Barrat, J., Feigelman, M., Kurchan, J. et at. (eds.) Slow 
Relaxations and Nonequilibrium Dynamics in Condensed 
Matter (Springer, 2004). 

11 Hedges, L. O., Jack, R. L., Garrahan, J. P. & Chandler, D. 
Dynamic order-disorder in atomistic models of structural 
glass formers. Science 323, 1309-1313 (2009). 

12 Garrahan, J. P. & Lesanovsky, I. Thermodynamics of 
quantum jump trajectories. Phys. Rev. Lett. 104, 160601 
(2010). 

13 Heyl, M., Polkovnikov, A. & Kehrein, S. Dynamical Quan- 
tum Phase Transitions in the Transverse Field Ising Model 
(2012). arXiv:1206.2505. 

14 Vogl, M., Schaller, G. & Brandes, T. Criticality 
in transport through the quantum Ising chain (2012). 
arXiv: 1208.5989. 



6 



15 Foss-Feig, M., Hazzard, K. R. A., Bollinger, J. J. & Rey, 
A. M. Non-equilibrium dynamics of Ising models with de- 
coherence: an exact solution (2012). arXiv: 1209.5795. 

16 Levitov, L., Lee, H. & Lesovik, G. Electron counting statis- 
tics and coherent states of electric current. J. Mat. Phys. 
37, 4845-4866 (1996). 

17 Esposito, M., Harbola, U. & Mukamel, S. Nonequilibrium 
fluctuations, fluctuation theorems, and counting statistics 
in quantum systems. Rev. Mod. Phys. 81, 1665 (2009). 

18 Lecomte, V., Appert-Rolland, C. & van Wijland, F. Ther- 
modynamic formalism for systems with markov dynamics. 
J. Stat. Phys. 127, 51-106 (2007). 

19 Flindt, C, C. Fricke et al. Universal oscillations in counting 
statistics. Proc. Natl. Acad. Sci. USA 106, 10119 (2009). 

20 Plenio, M. B. & Knight, P. L. The quantum-jump ap- 
proach to dissipative dynamics in quantum optics. Rev. 
Mod. Phys. 70, 101-144 (1998). 

21 Gardiner, C. & Zoller, P. Quantum noise (Springer, 2004). 

22 Barreiro, J. T. et al. An open-system quantum simulator 
with trapped ions. Nature 470, 486-491 (2011). 

23 Blatt, R. & Roos, C. F. Quantum simulations with trapped 
ions. Nature Physics 8, 277-284 (2012). 



Garrahan, J. P. et al. Dynamical first-order phase tran- 
sition in kinetically constrained models of glasses. Phys. 
Rev. Lett. 98, 195702 (2007). 

Ruelle, D. Thermodynamic formalism (Cambridge Univer- 
sity Press, 2004). 

Demboo, A. & Zeitouni, O. Large Deviation Techniques 
and Applications (Springer, 1998). 

Nazarov, Y. & Kindermann, M. Full counting statistics of 
a general quantum mechanical variable. Euro. Phys. J. B 
35, 413 (2003). 

Miiller, M., Hammerer, K., Zhou, Y. L., Roos, C. F. & 
Zoller, P. Simulating open quantum systems: from many- 
body interactions to stabilizer pumping. New J. Phys. 13, 
085007 (2011). 

Reimann, P. Foundation of statistical mechanics under 
experimentally realistic conditions. Phys. Rev. Lett. 101, 
190403 (2008). 

Bena, L, Droz, M. & Lipowski, A. Statistical mechanics 
of equilibrium and nonequilibrium phase transitions: The 
Yang-Lee formalism. Int. J. Mod. Phys. B 19, 4269 (2005). 



